THE IN-MEDIUM FEW-BODY PROBLEM 



S. A. SOFIANOSi AND M. BEYERi-^ 

^Physics Department, University of South Africa, 
P. O. Box 392, Pretoria 0003, South Africa 

Fachbereich Physik, Universitdt Rostock, 
18051 Rostock, Germany 



Abstract. We are concerned with few-particle correlations in a fermionic 
system at finite temperature and density. Within the many-body Green 
functions formalism the description of correlations is provided by the Dyson 
equation approach that leads to effective few-body equations. They contain 
the dominant medium effects, which are self energy corrections and the 
Pauli blocking. Hence the effective two-body interactions between quasi- 
particles are momentum/energy-dependent and therefore they can be us- 
esed in the medium modified, momentum space, integral AGS equations 
for three- and four-body systems. To investigate correlations and clusters 
beyond four-body, we employ, instead, the configuration space two- variable 
integro-differential equations (IDEA) for A-body bound systems which are 
based on Hyperspherical Harmonics and the Faddeev decomposition of the 
wave function in two-body amplitudes. This requires the transformation 
of the energy dependent two-body interactions to equivalent local, energy 
independent, ones. To achieve this we use inverse scattering techniques the 
resulting interactions being, on- and (to all practical purposes) off-shell 
equivalent to the energy dependent potentials. In this way we obtain bind- 
ing energy results for the 2-, 3-, 4-, and 16-particle in a medium at a 
finite temperature and various densities. Several aspects of the problem 
are discussed and the behavior of the potential surfaces obtained in the 
extreme adiabatic approximation, below and above the Mott transition, is 
investigated. 
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1. Introduction 

The properties of nuclear matter can be studied using various methods such 
as path integral methods, computer simulations, or perturbation theory 
using many-body Green functions [1]. In the latter, which will be followed 
here, the equations of state may be expressed in terms of many-body Green 
functions and the inclusion of correlations is provided by the Dyson equa- 
tion approach developed for this purpose in [2, 3, 4, 5, 6, 7], and reviewed 
in [8]. Within the Dyson equation approach the familiar in-medium two- 
body equations, known as Galitskii-Feynman or Bethe-Goldstone equations 
(depending on details), can be derived [1]. Recently, these equations have 
been extended to investigate correlations and clusters of more than two 
particles and applied to various aspects of fragmentation in a heavy ion 
collision at intermediate energies [9, 10, 11, 12, 13, 14, 15, 16, 17, 18]. For 
three- and four-body correlations this approach leads to modified momen- 
tum space integral equations of the Alt-Grassberger-Sandhas (AGS) [19, 20] 
type. 

For correlations with more particles (A > 4), relevant e.g. in fragmentation 
of heavy ion collisions at intermediate energies [21], the momentum space 
integral equations are not appropriate and thus one has to resort to effective 
two-body interactions for the clusters involved and consider the system 
under consideration as a few-body system. For instance, the ^^O may be 
considered as a A-a system interacting via the a — a interaction. In such a 
formalism, however, one faces the problem of how to construct the inter- 
cluster interactions which is a difficult task to implement. Apart from this, 
most of the existing and well established theories for the description of the 
medium-light nuclei, 5 < ^ < 20 say, are in configuration space and are 
based on energy-independent two-body potentials. Therefore, the energy 
dependent effective interactions can not be directly employed and thus one 
should resort to some localization procedure which transforms the V{E, r) 
for each chemical potential and temperature T, to an equivalent local 
potential (ELF) VL{r) which is energy- independent. The construction of 
the VL{r) is thus desirable for two reasons. Firstly, to investigate how the 
characteristics of the in-medium two-body interaction are manifested in the 
presence of the Fermi function f{E). Secondly, to enable us to employ a 
configuration space formalism describing the A-body system embedded in 
nuclear matter. 

In this respect, the two-variable integro-differential equation formalism 
(IDEA) in configuration space [22, 23, 24], valid for three- and many-body 

systems is employed. The IDEA is based on the Faddeev formalism, where 
the A-body wave fTinction can be written as a sum of two-body amplitudes, 
and on the Hyperspherical Harmonics. The reduced equations give bind- 
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ing energy results for the A-body which are in excehent agreement with 
variational, Green Function Monte Carlo, and other competing methods. 
In what follows, we shall describe first in Sect. 2 our formalism starting from 
a brief description of the effective in-medium equations, the localization 
procedure, and the IDEA method. Our results are discussed in Sect. 3 and 
our conclusions are summarized in Sect. 4 

2. Formalism 

Wc use the techniques of many-body Green functions. Reference to text- 
book treatments is given in [1]. The Dyson equation approach to cor- 
relations along with several approximation schemes has been developed 
in [2, 3, 4, 5, 6, 7], reviewed in [8] and will be briefly summarized in the 
next section. In recent years we have systematically applied this approach 
up to four-nucleon clusters [9, 10, 11, 12, 13, 14, 15, 16, 17, 18] and also 
generalized it to the light front to describe quark correlations in hot and 
dense relativistic nuclear matter [25, 26, 27, 28]. 

2.1. THE DYSON EQUATION APPROACH 

Let the many-particle Hamiltonian be given by 

H = Y, Ho{l, l')alav + ^ E ^2(12, l'2') a\ala2'ar (1) 

11' 121'2' 

where Ho is the kinetic energy, V2 a generic two-body potential, and a, 
are Fermionic destruction and creation operators. We define a chronological 
Green function for a given number of particles as follows 

gl/ = -i(TA„(t)4(t')) 

= -i {9{t - t'){Al{t')A^{t)) T 0{t' - t){Al{t')Aa{t))) . (2) 

The average (• • •) is taken over the exact ground state and the upper (lower) 
sign is for fermions (bosons) . The operators A^ {t) could be build out of any 
number of destruction and/or creation operators. For n-particle correlations 
Aa = 010203 . . .ttn- The time dependence of the operators is given in the 
Heisenberg picture by A{t) = e*^*Ae~*^*. Note that the multi-particle (or 
cluster) Green function is defined at equal or global time. 
In finite temperature formalism the above definition can be generalized. For 
a grand canonical ensemble, relevant here, the Heisenberg picture assumes 
the form A{t) = e^^~^^'>'^ Ae^^^~^^^^'^ , which includes the chemical poten- 
tial fi and the number operator N. For r = it (imaginary time) the original 
Heisenberg picture is formally recovered. The average is now taken over the 
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(equilibrium) grand canonical statistical operator p, viz. (• • •) = tr{p. . .}. 
To accommodate the convention of Fetter and Walecka [1] the imaginary 
time Green function is given by 

QV' = -{TrA^{r)Al{T')), (3) 

where the time ordering is according to the value of r, with the smallest at 
the right and the imaginary unit is dropped compared to (2). 
Dyson equations can be established for both forms [8, 10], 



. d 



The mass matrix that appears in (4) is given by 

Ml-/ = b{T-T')Ml^p^Ml:/ (5) 
{M^N)l^ = {\\A^.H\{t\A\{t)\) (6) 

{.MMJ-/ = E(^-[^",i?](r),[At,i7](r')])irreducible (7) 

7 

where we have used A/"^^ = ([A^, ^j3]±(T)). The first term in (5) is instanta- 
neous and related to the mean field approximation, the second term is the 
retardation term. This form suggests to first solve the mean field problem 
(neglecting the retardation term) and evaluate higher order contributions 
involving in the mean field basis. Hence (4) may be written in the 

following form 

^1:^0,^ = Sir-r')Kp + j:Ml^,g^,-/ (8) 

7 

or/ = QUp+nj jdfdf'gi:j,{N-'Mr%f gtr'- (9) 

7(5 

As a first approximation, we neglect the retardation part and hence solve 
only the mean field equation. 

It is convenient to use the Matsubara-Fourier representation [1] of the 
Dyson equations. This leads to the Matsubara frequency instead of the 
"time" derivative in (8) 

Att , , 

.A=-^+/X (10) 

where A = 2n for bosons and A = 2n + 1 for fcrmions. 

For a one body operator, (8) leads to the standard mean field approxi- 
mation (Hartree-Fock approximation) and an effective Hamiltonian that 
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describes uncorrelated quasiparticles. To be more specific, for A.,^ = ai the 
one-particle Green function, as follows from (8), the use of (6) and the of 
Matsubara-Fourier transformation, is 

G{z) = {z-e,r' (11) 

where the Matsubara frequency has been analytically continued zx z. 
The quasiparticle self energy is 

^-4+i:^=(12.r2,/. = ^ + AH^(l) (12, 

where 12 denotes proper antisymmetrization of particles 1 and 2. The Fermi 
function /j = /(ffj) for the i-th. particle in an uncorrelated medium is given 
by ^ 

(at a) = tr(pata) = /(e^) = ^^(,._^) ^ ^ ■ (13) 

The density operator of an uncorrelated medium is characterized only by 
single particle operators, viz. 

^° " tr(e-/^-^) ' ^ " ~ ^^^^ 

For a general operator Aa other than one body, (8) is called cluster mean 
field approximation. Using the density operator for an uncorrelated medium 
to evaluate the traces, the equation hierarchy decouples. This is a consistent 
way to arrive at equations for each cluster that can be solved with few-body 
techniques. 

The resolvent Go for n noninteracting quasiparticles is 

n 

Go{z) = {z- Ho)-^N = Ro{z)N, Ho = Y,ei (15) 

i=l 

where Gq, Hq^ and N are formally matrices in n particle space. The Pauli- 
blocking for n-particles is 

iV = /l/2.../n±/l/2.../„, 7=1-/ (16) 

Note that NRq = RqN. Straightforward but tedious evaluation of (8) using 
Aa = 010203 ... On (4) and the Wick theorem, leads also to the full resolvent 
G{z). For each number of particles in the cluster the resolvents have the 
same formal structure and after Matsubara-Fourier transformation can be 
written in a convenient way as in the isolated system, viz. 



G{z) = {z-Ho-V)-'N = R{z)N, V = ^^^2" ■ (17) 

pairs a 
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Note that V'^ ^ V and R{z)N / NR{z). To be specific, for an interaction 
in pair a = (12) the effective potential reads 

(123 . . . |iv(^'Vi'')|l'2'3' . . .) = (/1/2 - /i/2)F2(12, 1'2')533' ■ ■ ■ (18) 

For further use in the Alt-Grassberger-Sandhas (AGS) equations [19, 20], 
we give also the channel resolvent 

G^{z) = {z-Ho- N^V^r^N = Ra{z)N. (19) 

For the two-body case as well as for a two-body subsystem embedded in the 
n-body cluster the standard definition of the t matrix leads to the famihar 
t matrix equation for finite temperature and densities [1], 

T^iz) = + V^N^Ro{z)T^{z). (20) 

This equation, given in the quasiparticle approximation, is used as a driv- 
ing kernel for the few-particle equations. Going beyond the quasiparticle 
approximation would extend the scope of the present work. Note, however, 
that improvements can be systematically achieved by including the retar- 
dation part of the effective mass matrix, although technically this is quite 
a difficult task. If multi-particle correlations are neglected (binary collision 
approximation) it is possible to solve the t matrix in a self consistent way, 
see e.g. [29] . 

2.2. EFFECTIVE IN-MEDIUM EQUATIONS 

The in-medium Schrodinger type equation equivalent to (20) is given by 
[Ho - z)^{l2) + ^(1 - /i - 72)^2(12, 1'2')*(1'2') = 0. (21) 

1'2' 

In the present work we are interested in a region of rather low density in the 
vicinity of the Mott transition (~ /9o/10, where po = 0-16 fm~^ is normal 
nuclear matter density). We therefore use an effective mass approximation 
for the quasi-particle energy e; i.e. after evaluation of Eq. (12) we fit the 
effective mass rUeS via 

The constant shift A^^ can be absorbed in a redefinition of the chemical 

potential /Ucg = /i — A^^. To simplify our analysis further, we assume that 
the two-body system is at rest in the medium. It is well known that the 
influence of the medium fades away for larger relative momentum respective 
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to the medium [30]. Upon introducing relative and cm. coordinates for the 
two-particle system (p and P = 0), the kinetic energy term reads 

^ - 2 _ J- =-E (23) 



2meff 2meff 2meff 2 

and E = z — -Econt is the two-body binding/scattering energy in the cm. 
system. In the effective mass approximation and for P = 0, -Econt is given by 
-Econt = 2A[^^[30]. Using a local bare nuclcon nucleon potential V2{r), where 
r is the conjugate of p, the effective potential becomes energy dependent 

V{E,r) = {l-2f{E))V2{r), (24) 

where the Fermi function is now explicitly given by 

The effective in-medium Schrodinger-type equation for scattering states 
then reads ^ 

+V{E,r)-E\ il){E,r)={). (26) 



Note that this equation holds for E > Q, since E = p'^/rues for the eigen- 
values. 

The potential V{E, r) can be employed to study its effects and how is 
manifested in ^-body systems. For example, for A = 2, one may study its 
analytical behavior in the complex fc-plane as the nuclear density p changes 
[31] or study the influence it has on clusters embedded in the medium. 
Due to its energy-dependence, it can be employed in momentum space 
integral equations for three- [13] and four-body [16] systems embedded in 
the medium. 



2.3. ENERGY-INDEPENDENT INTERACTIONS 

Equivalence between two potentials Vi and V2, irrespective of their na- 
ture and structure, means that the scattering wave functions for the two- 
potentials are related by [32] 

ui{r) = f{r)u2{r) (27) 

with 

/(O,--! (28) 

i.e the two wave functions must be asymptotically the same and thus the 
phase shifts are also the same. In our case, however, we have to obtain 
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an equivalent potential to V{E, r) which should not only provide the same 
asymptotic scattering wave function but it should generate differences in 
the interior region which are minimal (/(r) ~ 1 V r) while at the same 
time the corresponding binding energy obtained by the ^^-independent 
interaction should be the same as the one obtained from V{E}),r). 
Prom all methods employed to construct ^'-independent ELP's, the inverse 
scattering technique is the most natural one. The two-body potential is, in 
this case, directly obtained from the available bound and scattering infor- 
mation in a unique way and without an a priori assumption about the 
shape and range of the interaction. This not only guarantees the recovery 
of the on-shell behavior of the wave function but of the bound states -when 
present- as well. The off-shell differences between the original wave func- 
tion and the one generated by the ELF potential obtained by inversion, 
described by the function /(r), Eq. (28), are also minimal. From all variant 
inverse scattering methods, the fixed-£ inversion of Marchenko [33, 34, 35] is 
perhaps the most suitable to employ in order to construct an ^^-independent 
ELF to V{E, r) as defined by (24). The required input is the S'-matrix Si{k), 
defined by 

Se{k) = = e^^^^W , (29) 

where fe{k) is the Jost-function the analytical properties of which describe 
the scattering, bound, and resonances of the system and <5^(A;) arc the scat- 
tering phase shifts at the energy E = h^k"^ /2fj, which are used to obtain the 
Fourier transform J^i(r,r') of the S-matrix Si{k), 

= ^J^h'^''\kr)[l-Se{k)]h^^\kr')dk 

+ J^A^eh^+\bnr)h'i^\bnr') , (30) 



n=l 



where h^^^ {z) is the Riccati-Hankel function, Nf, is the number of bound 

states, and An£ are the corresponding asymptotic normalization constants 
[35] . The J-'i in turns serve as an input to the Marchenko integral equation, 

/oo 
Ki{r,s)J'e{sy)ds = . (31) 

Once the K^^r, r') is obtained from the solution of (31), the potential V^(r) 
for the partial wave I is calculated from the relation 

W = -2^^. (32) 
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It is clear that while in the Schrodinger equation one starts from the knowl- 
edge of the potential to obtain the wave function which provide us with the 
physical information, in the inverse scattering procedure one uses instead 
experimental scattering data (or data obtained via a theoretical approach) 
to construct the underlying unknown potential first. The practical imple- 
mentation of the method has been described in [31, 36]. 

2.4. INTEGRO-DIFFERENTIAL EQUATIONS APPROACH 

The method we employ to study the dynamics of the A-boson-*^ like sys- 
tem is based on the Faddeev decomposition of the wave function and on 
Hyperspherical Harmonics. The A-body Schrodinger equation reads 



i<j<A 



*(x) = (33) 



where x stands for all the coordinates of all particles in the center-of mass 
frame, V{rij) is a central two-body potential, and rij = |xj — Xj|. Assuming 
that the wave function is written as a sum of amplitudes, 

^r(x)= E $(r,,-,x)=%^](x) Yl nrij,r) (34) 
i<j<A i<j<A 

where (x) is a harmonic polynomial of minimal degree and r is the 

hyperradius r = [(2/yl) J2i<j<A''^ijj obtains the Faddeev-type equation 
for the amplitude F(rij,r) 

[T - E]H^Lrr.]i^)F{rij,r) = -V{nj)H[L^]{^) ^ F{rki,r) (35) 

k<l<A 

This equation may be modified by introducing in both sides the average 
of the potential over the unit hypersphere i.e the so-called hypercentral 
potential Vb(r) defined by [22] 

Voir) = V{r^{l + z)/2)W{z)dz (36) 

where z = 2rfj/r^ — 1 and where the weight function W{z) is given by 

W{z) = {l~zr{l + z)f^ (37) 

^Presently we neglect spin-isospin dependence, which means that the effective few- 
body equations derived are for bose-Uke particles. This restriction seems not severe as it 
might be interpreted as a spin averaging. This kind of averaging procedures (e.g. angular 
momentum averaging of Pauli factors) are widely used for quasipaxticles and lead for 
example to the possibility of angular momentum expansion. This does not affect the 
Fermi statistics in use. 
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with a = (L> - 5)/2, /3 = 1/2, and D = 3{A-1). Then instead of (35) one 
has the modified Faddeev equation 

l^(ri,)-^o(r)j%^](x) J2 <^(r«'0 (38) 
-' k<l<A 

It is noted that by summing aver all pairs (ij) one recovers the Schrodinger 
equation. Letting (/)(rjj,r) = P{z,r)/A^~^^^'^ for the modified Faddeev 
component, multiplying from left by H^^^ ](x), and integrating over the 
(N — 1) -hyper angles dOjv-i (except those of the (z,j)-pair) we obtain the 
integro-differential equation (IDEA) for the amplitude P{z, r) [22, 23] 





I 4 1 


{- 




t m 





+ 



^(^-1) 



Vo{r)-E^ P{z,r) 

V(rsJ(l + z)/2) - Voir) 



n(z,r) (39) 



where we use the abbreviations 
52 C{C + 1) 



and 



n(z,r) = P(z,r)+ / f{z,z')P{z',r)dz' . 



In the above, £ = (D — 3)/2 and the function f{z, z') is given by 

f{z, Z') = W{Z') Y,Uk - l)PK{z)PK{z')/hK 
K 



(40) 
(41) 

(42) 



where Pk{z) = P^'^{z) are Jacobi polynomials associated to the weight 
function W{z). The constant fx is given by 



1 



2iA - 2)Pk{-1/2) + (-4 2)(^ ^^ P^(-l) 



while /ix, a normalization constant, by 



hK = j ^ [PK{z)fW{z)dz 



IPk{1) (43) 



(44) 



Details on the IDEA formalism can be found in Refs. [22, 23] and will not 
be discussed here. 
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In the adiabatic approximation one assumes that the radial motion (coor- 
dinate r) and the orbital motion (coordinate z) are nearly decoupled i. e we 
may right P{z,r) ~ Px{z, r)ux(r) and that P\{z,r) is varied slowly with r 
and thus one can omit its derivatives with respect to r. Then Eq. (39) can 
be split into two equations {h^ /m = 1) 



- V>, + Ux{r) 



Px{z,r) 



V{r^{l + z)/2) - Voir) 



Uxiz,r) (45) 



and 



Dr + 



AiA - 1) 



Voir) + Uxir)) 



^caa 



u\[r, 



= 



(46) 



The U\ir) is the so-called eigen-potential sometimes referred to as potential 
surface. It is a characteristic for each system and it can be used to calculate 
reliably the yl-boson like bound states via (46). The solution thus obtained 
provide us a lower bound for the binding energy Ef^'^ [22, 23, 24]. 

3. Results 

In our investigations we employ the Volkov potential [37] which is widely 
used in model nuclear structure calculations, 

V{r) = -83.3400196 cxp(-(r/1.6)^) + 144.843409 cxp(-(r/0.82)2) (47) 

the two-body binding energy being E2 = 0.5462 MeV. The medium-modified 
energy dependent potentials ViE,r) are constructed using Eq. (24). Em- 
ploying ViE, r) one can obtain phase shifts at any energy E corresponding 
to a specific and effective mass rrics. In the present work we use two 
nuclear densities namely p = 0.003 fm""^ and p = 0.034 fm~^ and T = 10. 
The low energy behavior of these phase shifts are shown in Fig. 1. The 
distortion of the potential for p = 0.003 fm~^ results in different phase 
shifts but the potential still sustains a bound state at E2 = 0.0282 MeV. 
In contrast the phase shift for p = 0.034 fm~^ start from zero and thus, 
according to Levinson's theorem, no bound exists in this case. 



"<5 



180- 


1 




120 






60 






n 


1 



2.5 
E [MeV] 



Figure 1. The low energy behav- 
ior of the phase shifts without dis- 
tortion (— ), for p = 0.003 fm"^ 
(---), and for p = 0.034 fm-^ (- •-) 
at T = 10 MeV. The latter phase 
shifts start from zero implying the 
absence of a bound state. 
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T 




5 10 5 10 

r [fm] r [fm] 

Figure 2. Comparison of the scattering wave functions V"(r) at E = 25McV obtained 
with the ^-independent Marchenko potential Vi,{r) ( — ) with those of V{E,r) ( — ) for 
p = 0.003 &n-3 (Fig. a) ) and p = 0.034 frn'^ (Fig. b) ) at T = 10 MeV. 

Using the above phase shifts we constructed the ^^-independent potentials 
by means of the Marchenko inverse scattering method, the quahty of which 
can be seen in Fig. 2 where the original and the resulted wave functions at 
E = 25 MeV are drawn. It is seen that, to all practical purposes, the two 
wave functions are the same i.e the potential V{E,r) is not only on-shell 
equivalent but (nearly) off-shell equivalent to the one obtained by inversion, 
VL{r), as well. Similar results are obtained at all energies. 
Once the quality of the interactions is checked, they are employed in the 
IDEA equations (45) and (46) to obtain the binding energies for the three-, 
four- and sixteen-body systems in the adiabatic approximation. The bind- 
ing energy results for the three densities used are given in table 1. It is seen 
that with the bare interaction used, the three-body is the first to 'dissolve' 
in the medium for density p = 0.034 fm^"'. The four-body results show sim- 
ilar trends to the three-body system. However, the binding energy for the 
16-body system, is not as drastically reduced and for p = 0.034 fm~^ which 
is quite higher than the Mott density, the binding energy is still significant. 
In Fig. 3 the ELP's for the medium modified 2-body interaction is plotted 



TABLE 1. Binding energy results for 2-, 3-, 4-, and 16-body 
systems. 



p (fm-3) 


E2 (MeV) 


Ea (MeV) 


E4 (MeV) 


E16 (MeV) 





0.5462 


8.6704 


30.7250 


1630.1462 


0.003 


0.0282 


4.7202 


22.2798 


1554.0695 


0.034 






5.4093 


1410.5321 



together with the potential surfaces obtained from the solution of (45). It is 
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seen that the essential features of the medium modified two-body interac- 
tion with the Volkov (bare) potential are similar to those obtained Ref. [31] 
with the Malfliet-Tjon potential: As the density of the medium increases 
the attractive well becomes narrower and eventually for densities above the 
Mott one, the two nucleons become unbound. It is interesting to note that 
at these densities the repulsion in the interaction region appears. The ana- 
lytical properties for the corresponding Jost function of the potentials have 
been studied in Ref. [31] where it was shown that the bound state, instead 
of being continued off the imaginary axis to the resonances region, moves 
along the negative imaginary fc-axis, i.e. it becomes an anti-bound state - 
at least for the densities used where the repulsive hump in the interaction 
region is weak. The potential surfaces exhibit similar behavior and in ad- 
dition they become shallower. The eigen-potentials for the 16-body system 
tend, in addition, to shift the repulsion in the interior region in a more 
pronounced manner. 




Figure 3. The two-body potentials and the corresponding potential energy surfaces (in 
the extreme adiabatic approximation) for 3-, 4- and 16-body systems for zero distortion 
(- - -), for p = 0.003 fm-3 (---), and p = 0.034 fm"^ (-•-). 
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4. Conclusions 

We propose a method to investigate ttie A-body system embedded in the 
medium at a finite temperature and for densities below and above the Mott 
density. Unlike previous investigations which were carried out in momentum 
space, the present one is in configuration space and thus advantageous in 
visualizing the effects of the medium on the system under consideration. 
To achieve this, we firstly apply inverse scattering techniques to transform 
the medium modified two-body interactions which are energy dependent, 
into an energy independent potentials. This allow us to employ the potential 
directly in configuration space where various formalisms exist. One such 
formalism, is the IDEA method for few- and many-body systems which, 
apart from its simplicity, gives results which are in excellent agreement 
with those of other methods and therefore they can be safely employed in 
our investigations on the in-medium few-body problem. 
We demonstrate our approach by calculating the binding energies for the 
2-, 3~, 4-, and 16-boson like systems at finite temperatures and densities. 
The behavior of the adiabatic approximation potential surfaces, below and 
above the Mott density, is also investigated. We found that the adiabatic 
potentials exhibit similar features to those found in Ref. [31], namely, that 
the attractive well moves towards the interior region while at the same time 
is becomes narrower as the densities become larger. The analytical prop- 
erties of the corresponding Jost function in the complex /c-planc for these 
potentials were investigated in [31] where it was shown that the bound state 
moves, downwards on the imaginary fc-axis and eventually 
become an anti-bound state instead of a resonance (as implied by the ap- 
pearance of a small repulsive hump in the interaction region). This reflects 
a vanishing imaginary part of the respective two-body spectral function. 
Hence, the two-body cluster (even for 6 < 0) retains the character of a quasi- 
particle, i.e. the correlation does not decay. In a more general approach the 
self-consistency requirement [29] as well as three-body collisions [11] lead 
to imaginary parts of the two-body spectral function, and therefore to a 
finite width (life time) of the two-particle cluster irrespective of b smaller 
or larger than zero. 

It is interesting to note that the change of the eigen-potential for the 16- 
body shows a lesser sensitivity to the density changes and the binding 
energies are still quite deep. Further investigations with other bosonic sys- 
tems and two-body forces are needed to elucidate the behavior and phase 
transition of large clusters in the medium. 
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